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ABSTRACT 

We investigate the impact of gas evolution on the expected yields from deep Sunyaev-Zel'dovich (SZ) 
effect surveys as well as on the expected SZ effect contribution to fine scale anisotropy in the Cosmic 
Microwave Background. The approximate yields from SZ effect surveys are remarkably insensitive to gas 
evolution, even though the observable properties of the resulting clusters can be markedly different. The 
CMB angular power spectrum at high multipoles (£ > 2000, < 5') due to the SZ effect from clusters is 
quite sensitive to gas evolution. We show that moderate resolution (< 1') SZ effect imaging of clusters 
found in deep SZ effect surveys should allow a good understanding of gas evolution in galaxy clusters, 
independent of the details of the nature of the gas evolution. Such an understanding will be necessary 
before precise cosmological constraints can be set from yields of large cluster surveys. 

Subject headings: cosmic microwave background — cosmology: theory 
universe — cosmological parameters 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect has recently become 
a valuable observational tool (Birkinshaw 1999). In the 
past few years, observational programs have shifted from 
attempting to detect the effect to using the SZ effect to 
map out massive clusters of galaxies, study the intra- 
cluster medium (ICM) , and constrain cosmological param- 
eters (Carlstrom et al. 2000). 

The thermal SZ effect is a distortion of the cosmic mi- 
crowave background (CMB) spectrum caused by hot gas 
along the line of sight to the surface of last scattering. 
The cool CMB photons undergo Compton scattering on 
the hot electrons, gaining on average a small amount of 
energy in the process, creating an intensity decrement at 
low frequencies (y < 218GHz) and an increment at high 
frequencies. 

Current instruments are now regularly detecting and 
imaging clusters at high signal-to-noise, and the next gen- 
eration of instruments should be capable of mapping fairly 
large portions of the sky as a means of finding clusters of 
galaxies. Several works (e.g., Bartlett 2000; Holder et al. 
2000; Barbosa et al. 1996; Kneissl et al. 2001) have pre- 
dicted the number of clusters that could be expected in 
future SZ surveys and their angular power spectrum. The 
survey yields are quite impressive; the next generation of 
SZ instruments should be able to detect several clusters 
per day. 

In this work, we take a slightly different perspective. 
Rather than asking how many clusters an SZ survey will 
find, we instead focus on which clusters we would like to 
find and what we can learn from these objects. We use 
simple models for the gas distribution of galaxy clusters 
to investigate the impact of non-gravitational heating on 
observable cluster properties, focusing on the SZ effect. 
Much work has gone into using the X-ray emission from 
nearby groups and low-mass clusters to constrain feed- 
back from galaxy formation and non-gravitational heat- 
ing (e.g., Bower et al. 2000; Cavaliere, Menci, and Tozzi 



1999; Balogh, Babul, and Patton 1999) or alternatively the 
preferred cooling of low-entropy gas (Bryan 2000). The 
effects of feedback or cooling should be even more impor- 
tant at higher redshift and for low masses, where the SZ 
effect will be an important probe. 

In §2 we outline our models and the important SZ ob- 
servables. We show the expected SZ properties of clusters 
in §3 and angular power spectra for these models in §4. 
Finally, in §5, we discuss our results and implications for 
planned SZ surveys. 

2. SZ EFFECT FROM GALAXY CLUSTERS 

X-ray observations indicate that the cluster gas is ap- 
proximately isothermal (Irwin and Bregman 2000), and 
that the gas density profile can be approximated by a 
"beta model" (Mohr, Mathiesen, and Evrard 1999): 

where p gfi is the central gas density and r c is a core radius. 
We choose to fix (3 = 2/3, with the knowledge that this 
will give a reasonable SZ profile even if the true SZ best- 
fit value of (3 for clusters is slightly higher (Reese et al. 
2000; Grego et al. 2001). With this choice of f3, the mass 
is divergent, and must be truncated at some radius. We 
set the truncation radius to be the virial radius, derived 
from the spherical collapse model, as outlined below. 
For this cluster model, the central decrement of the clus- 
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where y is the Compton y parameter, R v is the virial 
radius, T e is the electron temperature, n is the central 
electron number density, <tt is the Thomson cross-section, 
and f{v) is a dimensionless frequency factor (e.g., see 
Birkinshaw 1999 ), which we will take to be -2 (i.e., the 
Rayleigh- Jeans limit). The central decrement is redshift- 
independent for a given set of cluster properties, allowing 
a distance-independent probe of cluster properties. 
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An observable which is more readily interpreted is the 
integrated SZ flux density from clusters. Assuming that 
the cluster is approximately isothermal and noting that an 
integral over solid angle can be written as dSl = dA/d A , 
we can write 



s = s i ;>_^Z_ = Sof{v)^a T A ' 
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where N e is the total number of electrons, dA is the an- 
gular diameter distance and S Q gives the conversion from 
temperature to flux density. At a frequency of 30 GHz, the 
conversion between flux and Y = J ydfl is approximately 
S = 12. 7 (Y/ar cm 2 ) Jy. The mass distribution does not 
matter in this case, as long as the integration over solid 
angle extends to the cluster virial radius and the cluster is 
roughly isothermal. 

The simple relation between the total flux density and 
the total mass and electron temperature make the SZ a 
particularly valuable tool for searching for clusters. Clus- 
ter detections are largely decoupled from the details of 
the structure of clusters, allowing an unbiased sample for 
studies of both cosmology and cluster properties at any 
redshift. 

2.1. Mass-Temperature Relation 

In connecting with SZ observables, a relation is needed 
between mass and temperature. Such a relation can be 
derived from the spherical collapse model (Lahav et al. 
1991) or can be found from numerical simulations. We 
have chosen the hybrid approach of Bryan and Norman 
(1998), where the functional form is taken from the spher- 
ical collapse model and the normalization is found from 
simulations. 

Using a notation where H(z) = 100hE(z) km s _1 Mpc -1 , 
we use 

T gas = 1.4 x 10 7 (A(z)£(z) 2 )sa4 K , (4) 

with Mi 5 the virial mass in units of 1O 15 /i _1 M , A(z) is 
the mean cluster density relative to the critical density at 
that redshift and the virial radius is defined as the radius 
that encloses this mean density. 

2.2. Evolution of Cluster Gas 

The SZ observables will be sensitive to the density pro- 
file of the ICM, so it is important to model the ICM dis- 
tribution both as a function of mass and as a function 
of redshift. The simplest evolution is self-similar evolu- 
tion, where the central density scales with the background 
density and sizes will be proportional to the cosmic scale 
factor. In the absence of physics other than gravity, we 
would expect something like this sort of evolution, since 
there is no preferred scale. 

It is wcll-known(e.g., Kaiser 1991; Evrard and Henry 
1991; Ponman, Cannon and Navarro 1999) that X-ray ob- 
servations suggest that gas evolution is not self-similar. If 
reionization added significant entropy to the gas or if there 
has been substantial feedback into the ICM from galaxy 
formation, the excess entropy would shift the ICM to a 
higher adiabat, leading to a final state of the cluster that 
is less compact. In such a situation, the outer parts of 
the cluster might be expected to be largely unchanged, 
but the details of the cluster profile might depend on the 



details of the entropy injection (or preferential removal of 
low entropy gas). 

We will assume that the outer parts of clusters are un- 
affected by feedback or preheating, and approximate the 
dark matter density profile as a singular isothermal sphere, 
while the gas density has a well-defined core that is set by 
the minimum entropy of the gas. For self-similar evolution, 
this central minimum entropy is set by gravitational pro- 
cesses, while preheated models will have extra entropy in- 
jected by other processes, which may or may not be larger 
than the entropy from gravitational heating. 

As a simple model, we adopt the notion of an entropy 
floor (Ponman, Cannon, and Navarro 1999), where feed- 
back has provided a uniform amount of entropy, and grav- 
itational heating has provided an amount of entropy that 
increases with increasing mass. This is a slight modifi- 
cation of the constant entropy core model of Evrard and 
Henry (1991). With an entropy floor rather than a fixed 
central entropy, high-mass clusters can have a central en- 
tropy which is significantly higher than the entropy floor, 
more in line both with physical intuition and observations. 

We define a measure of the entropy as 
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and take the central entropy to be the entropy expected 
from self-similar evolution plus a constant entropy floor. 

The effect of the central entropy floor is to set a max- 
imum central density for a given temperature. For high- 
mass clusters, the central entropy is set by the gravita- 
tional heating; for low-mass clusters, the entropy from 
gravitational heating is relatively small and the central 
entropy is determined by the entropy floor, leading to a 
temperature-dependent central density. 

2.2.1. Self-Similar Evolution 

For self-similar evolution (sfi oor = 0), the core radius 
should be a fixed fraction of the virial radius. We choose 
this ratio 1Z = R v /r c = 10, in rough agreement with obser- 
vations of clusters (Mohr, Mathiescn, and Evrard 1999). 
Under the assumption that the gas density at the virial ra- 
dius is equal to the global baryon fraction, fs = f2_B/T2 m , 
times the mass density at the virial radius, we solve for 
the central density: 
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In general, the global baryon fraction is not equal to 
the cluster gas mass fraction. Simulations without feed- 
back from galaxy formation typically find values for the 
cluster gas mass fraction that are only slightly lower than 
the input global baryon fraction (Evrard 1997). As we 
outline below, it is possible for the cluster gas mass frac- 
tion to deviate strongly from the global baryon fraction if 
non-gravitational heating is introduced. 

The entropy expected from self-similarity for a cluster 
of a given mass is easily obtained from equation 5 using 
equations 4 and 6, and can be seen to simply scale with 
the virial temperature. 

The central SZ decrement assuming self-similar evolu- 
tion and solar abundances of hydrogen and helium (ji e = 
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1.14), from equation 2, is 

AT = 295 [iK f B h arctanft Mi 5 Ai 7 $E{z) 2 

(7) 

The total flux density, assuming an observing frequency of 
30 GHz, can be calculated from equation 3: 

S = 146 mJy f B h M^ 3 (A 178 E{zf) 1/3 
(I + 1Z 2 )(1Z- tan' 1 K) / d A {z) y* 

where A178 = A/178. In a universe with Q m = 1 the 
spherical collapse model predicts A = 178, with smaller 
values for low-density universes. 

2.2.2. "Preheated" Models 

For a given temperature, the central density is fixed by 
the central entropy, and we can use equation 6 to solve 
for the appropriate value of the core-to- virial ratio TZ. For 
the self-similar case the central entropy simply scales with 
T and TZ = constant. A non-zero entropy floor breaks 
the self-similarity and TZ becomes a function of tempera- 
ture. We assume that the entropy from gravitational heat- 
ing is equal to that expected in the self-similar case and 
add entropy from non-gravitational heating. We can still 
use equations 7 and 8, but we must now solve for 1Z as 
a function of mass. For low-mass clusters, this leads to 
significantly different evolution and appearance, as well as 
a depressed cluster gas mass fraction. As an example, for 
= 1, the self-similar case would predict that the central 
decrement would scale &s M(l + z) 3 , whereas the entropy 
floor would predict a scaling as M 3 / 2 (l + z) 9 / 4 . This re- 
lation is different in both its mass scaling and its redshift 
dependence. In this scenario, low-mass high-redshift clus- 
ters will be significantly less compact than in a self-similar 
picture. This would have observable consequences, both 
in the properties of high-redshift clusters and in their sig- 
nature in CMB experiments (see §4). 

For clusters well below ~ 10 14 /i _1 M© the core radius 
can become significantly larger than the virial radius for 
high values of the entropy floor. In such cases, the gas 
distribution is still truncated at the virial radius in our 
models. This is simply indicating that the entropy of the 
central gas is sufficiently high that the central density is 
not significantly higher than the gas density near the virial 
radius. In this regime, the accretion process could be seri- 
ously affected by the entropy floor and it is not clear that 
our simple model is applicable. 

In this work, we assume a relatively high value for the 
entropy floor of s fi oor = 200 keV cm 2 as an extreme case. 
Current estimates of the entropy floor (Ponman, Cannon, 
and Navarro 1999) are roughly Sfi oor = 100 keV cm 2 . 

2.3. Cosmological Evolution of Cluster Abundances 

We modeled the cosmological evolution of cluster abun- 
dances with the Press-Schechter prescription (Press and 
Schechter 1974; Bond et al. 1991), which gives the comov- 
ing number density as a function of both mass and redshift. 
We followed the procedure outlined in Holder et al. (2000), 
with two minor modifications. The power spectrum was 
computed using the fitting functions of Eisenstein and Hu 
(1999), and the redshift evolution of the power spectrum 



was evaluated numerically from linear theory (e.g., Peebles 
1980). 

Given the comoving number density, it is straightfor- 
ward to obtain the number of clusters per steradian above 
some mass threshold Mu m as a simple integral over the 
comoving number density. To estimate the angular power 
spectrum of the thermal SZ, we follow Cole and Kaiser 
(1988), and assume that only Poisson contributions to the 
angular power spectrum are important. 

For this work, we assume that Q m — 0.3, Q\ = 
0.7, as — I and h = 0.65. The expected counts are 
very sensitive to cosmological parameters, as many au- 
thors have found. 

3. PROBING HIGH-REDSHIFT CLUSTERS 

The effects of preheating on cluster structure can be 
seen in Figure 1. The largest effects of preheating are 
seen at low masses and high redshift. These models have 
been designed to agree with observed properties of X-ray 
emitting clusters and groups at low redshift, so the true 
test of these models will be how they fare at high red- 
shift. In particular, important information can be learned 
about cluster-to-cluster variations in the amount of non- 
gravitational heating and also in the redshift evolution of 
the amount of preheating in clusters. 
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Fig. 1. — Central decrement and core radius for several masses at 
several redshifts (z = 0.5, 1, 1.5,2) for a preheated model (s/; oor = 
200keV cm 2 ) and self-similar evolution (s fi oor = 0). Mass is marked 
by open squares along each line and increases along each line from 
bottom to top in factors of 2, from 0.5 — 64 X 10 14 /i — 1 Mq. Lines of 
constant redshift are marked with similar symbols. Solid (open) tri- 
angles mark the limiting mass at which dN/dz > 1(0.01) per square 
degree per unit redshift for a ACDM cosmology. 

In Figure 1, we have marked the limiting survey mass 
at each redshift which is required to obtain a surface den- 
sity per redshift bin of either 1 or 0.01 clusters per square 
degree per unit redshift. In order to get a large sample 
of clusters at high redshift, the survey mass limit must be 
relatively low. 

Typical expected mass limits for SZ surveys are shown 
in Figure 2. PLANCK is expected to be able to survey 
most of the sky, at a resolution of ~ 5' down to a level 
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of ~ 10[iK. Deep ground-based surveys should be able to 
reach a comparable temperature uncertainty on a few tens 
of square degrees at a resolution of ~ 1'. 

The SZ is well suited for surveying, with the nice fea- 
ture that the survey mass limit should be largely decoupled 
from issues of preheating, as illustrated in Figure 2 by the 
small difference in the mass limits between the self-similar 
model and the fairly extreme preheated model. 




Fig. 2. — Differential number coants of clusters per square de- 
gree and mass limits (inset) for SZ surveys for no preheating (solid 
curve) and significant preheating (dashed curve) for a survey similar 
to PLANCK (bottom set of lines for dN/dz and top for mass limits) 
and deep ground-based surveys (top for dN/dz and bottom for mass 
limits). Flux density units are at 30 GHz (S ~ 12. 7 (Y/ar cm 2 )). 
Solid curves assume the ACDM cosmology outlined in the text, 
whereas the dotted curves show self-similar models with Q m = 0.33 
and erg = 0.9 and 0.95 for the top and bottom lines in the main 
panel, respectively. 

Rcdshifts for the clusters will be important to learn 
about ICM preheating, but it is interesting that a min- 
imum core radius is imposed by preheating which is not 
very sensitive to redshift. The simple existence of an ap- 
proximate minimum core radius would be a "smoking gun" 
for preheating, and the value of this core size is an indi- 
cator of the approximate value of the entropy floor. For 
example, using Sfi oor — 100 keV cm 2 will result in a min- 
imum core radius of roughly 30". Redshift evolution of 
the amount of preheating should be observable in the SZ 
properties of the cluster catalog. 

Of particular importance is that preheating should not 
significantly affect the number of clusters discovered in an 
SZ survey, only the appearance of the discovered clusters. 
However, as indicated by the dotted line in Figure 2, the 
change in expected yields is comparable to the change in 
yields expected from a shift in Q m by 10%. Clearly, non- 
gravitational heating will be an important consideration 
when trying to make precise estimates of cosmological pa- 
rameters from SZ surveys (e.g., Haiman, Mohr, and Holder 
2001). 

4. SMALL SCALE CMB ANISOTROPY 

The structure of high-redshift clusters can have impor- 
tant implications for studies of CMB anisotropies. At 



small angular scales, the angular power spectrum can be- 
come dominated by the thermal SZ effect. The magni- 
tude of the thermal SZ power spectrum is most affected 
by preheating at exactly the angular scales where it be- 
comes dominant, as shown in Figure 3, as also found in 
other studies (Holder and Carlstrom 1999; Komatsu and 
Kitayama 1999; Springel, White, and Hernquist 2000). 
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Fig. 3. — Angular power spectru^n from SZ clusters. Top solid 
curve shows the power spectrum due to all SZ clusters in the self- 
similar model, while the top dashed curve shows the power spectrum 
from clusters assuming a model with significant preheating. The 
middle solid curve shows the residual power spectrum for the self- 
similar model after all SZ clusters brighter than 15 mJy (at 30 GHz) 
have been removed, with the middle dashed line showing the pre- 
heated model. The bottom solid line and dashed curves correspond 
to the self-similar and preheated models, respectively, after subtrac- 
tion of all SZ clusters with flux brighter than 0.5 mJy. The expected 
contribution from primary anisotropies (unlcnsed) is shown as the 
light dotted curve, generated by CMBFast (Seljak and Zaldarriaga 
1996). 

Subtracting out the bright clusters can significantly re- 
duce the power at low multipolcs, but does not signifi- 
cantly affect the angular power spectrum at high I, where 
the thermal SZ effect is expected to exceed the primary 
anisotropies. The power at high i is mainly due to low- 
mass clusters at higher z, where non-gravitational heating 
can be very significant. The higher densities at higher red- 
shift lead to a lower self-similar prediction for the central 
entropy, making clusters at higher redshift especially sen- 
sitive to excess entropy. Deep SZ surveys will allow studies 
of the clusters which are contributing the majority of the 
signal at small angular scales. Shallow surveys will detect 
the clusters which provide the bulk of the power at larger 
angular scales, where the primary anisotropies are more 
than an order of magnitude larger than the thermal SZ 
effect, but will have very little information on the small 
angular scales at t > 2000. 

5. DISCUSSION AND CONCLUSIONS 

We have shown that the SZ effect is well-suited for stud- 
ies of non-gravitational heating. A well-planned survey 
should not have a selection function which strongly de- 
pends on the amount of non-gravitational heating, while 
SZ images of the resulting catalog will be very sensitive 
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to this energy injection. As a rcdshift-indepcndent probe, 
the SZ is particularly well-suited for high-redshift work, 
where non-gravitational heating can be very important. 

We have shown that the expected yield and redshift dis- 
tribution of clusters for an SZ survey is only weakly sensi- 
tive to preheating, but that this sensitivity is significant if 
one wishes to do "precision cosmology" with SZ surveys. 
This is true for either deep or shallow surveys. While high- 
mass clusters are less sensitive to preheating, the survey 
yield is more sensitive to the mass limit at higher masses. 
Thus we have a situation where it will be difficult to use 
high-mass clusters to learn about preheating, but we must 
understand it to do cosmology. For low-mass clusters, the 
mass limit is still important for the survey yield, but direct 
imaging of the resulting catalog of clusters should yield a 
wealth of information on preheating. 

The peak of the SZ angular power spectrum is very sen- 
sitive to the details of non-gravitational heating and direct 
images of high redshift, low mass clusters may be the best 
way of understanding this signal. Approximately half of 
the power at I > 2000, where the thermal SZ signal is 



stronger than the primary anisotropies, should be coming 
from clusters that could be detected in deep SZ surveys. 

Upcoming CMB experiments such as MAP and 
PLANCK, with their large beam sizes, arc not well- 
suited for finding high redshift and low mass clusters. 
Ground-based bolometer arrays and interferometers with 
sub-arcminute resolution will be well equipped for find- 
ing and studying these clusters. Deep SZ surveys will 
happen within the next few years and the information 
that they yield will be very valuable for constraining non- 
gravitational heating through galaxy formation or reion- 
ization. This information should be remarkably unbiased, 
and improvements in SZ imaging over the next several 
years should continue to improve our understanding of 
these processes. 

This is work supported by NASA LTSA grant num- 
ber NAG5-7986 and the DOE at Chicago. We acknowl- 
edge many useful discussions with Joe Mohr and thank U. 
Seljak and M. Zaldarriaga for making CMBFast publicly 
available. 
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